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We calculate the high energy neutrino spectrum from gamma-ray bursts where the emission arises 
in a dissipative jet photosphere determined by either baryonically or magnetically dominated dy- 
namics, and compare these neutrino spectra to those obtained in conventional internal shock models. 
We also calculate the diffuse neutrino spectra based on these models, which appear compatible with 
the current IceCube 40+59 constraints. While a re-analysis based on the models discussed here 
and the data from the full array would be needed, it appears that only those models with the most 
extreme parameters are close to being constrained at present. A multi-year operation of the full 
IceCube and perhaps a next generation of large volume neutrino detectors may be required in order 
to distinguish between the various models discussed. 



I. INTRODUCTION 

Gamma-ray bursts (GRBs) are a potential source of 
astrophysical high energy neutrinos which are currently 
being investigated with IceCube. The standard GRB in- 
ternal shock scenario of neutrino production [TJ [5] used so 
far to compare against the IceCube 40 string and 40+59 
string observations [3 , 4_ assumed some simplifications in 
the neutrino physics. Also, the astrophysical model itself 
of the GRB prompt gamma-ray emission based on the 
same internal shocks has been the subject of discussions 
in the gamma-ray community [SHE] , due to issues with 
the radiation efficiency and the spectral properties in the 
standard version of this internal shock scenario. 

For this reason, modified internal shock models that 
address these issues (e.g. |9Tll2]) as well as alternative 
models where the prompt gamma-ray emission arises in 
the jet photosphere have been considered (e.g. |13H18| . 
In such photospheric models the high radiative efficiency 
is due to dissipation processes in it. 

A separate question that has also been the subject of 
debate in the astrophysics community is whether the jets 
in such relativistic sources are dominated by the baryons 
or by magnetic fields, which imply different macroscopic 
acceleration rates, different proper densities in the jet 
rest-frame, and implying a major role for magnetic dissi- 
pation in the process of particle acceleration. Such mag- 
netically dominated jets in GRBs have been investigated 
by |H)| 22 and others, and the gamma-ray emission is 
ascribed in such models, again, to dissipative processes 
mainly in the photosphere, e.g. |23[ [24] . 

It is unclear at present whether the above mentioned 
modified internal shocks or the photospheric models are 



best for interpreting the prompt gamma-ray emission, 
nor whether the jet dynamics is dominated by baryonic 
or magnetic stresses (e.g. |25]). However, the expected 
neutrino emission is strongly dependent on the specific 
overarching dissipation and dynamic model of GRBs. For 
this reason, here we explore the neutrino features of the 
three main types of models which are currently under 
consideration. We illustrate the variety of astrophysical 
uncertainties involved in these models, and how these 
can affect the expectations for detection with IceCube or 
future instruments. 



II. THE DISSIPATIVE PHOTOSPHERIC 
SCENARIOS 

In the typical GRB model a high energy-to-mass ratio, 
jet-like relativistic outflow is launched which initially ac- 
celerates with a bulk Lorentz factor T, averaged over the 
jet cross section, whose dependence on distance from the 
center of the explosion can be parametrized as 

r = (r/r )". (1) 

This behavior is assumed valid up to a saturation ra- 
dius r sat = rQij 1 ^ 1 , where the Lorentz factor has reached 
the asymptotic value r sat = 77, where r/ ~ L/Mc 2 is the 
dimensionless entropy of the outflow, L and M being the 
average energy and mass flux. The index 1/3 < fx < 1 
ranges between the extreme \x — 1/3 magnetically dom- 
inated radial outflow and the usual /i = 1 baryonically 
dominated outflow regimes, e.g. |26| . In the extreme 
magnetic case and baryonic cases, the saturation radius 
is given by 
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(2) 

In the dissipative photospheric scenario, a fraction of 
the outflow bulk kinetic energy is converted into radiation 
energy via some dissipation mechanism in the neighbor- 
hood of the photospher^j] , giving rise to the "prompt" 
photon luminosity L 1 — e e L tot , where L tot — 10 53 L53 
erg/s is the isotropic equivalent total luminosity of the 
jet[j The photospheric radius is estimated by setting the 
Thomson optical depth r 7e ~ n. e oxilph/r = L where 
n' e = n' p L tot / An R 2 ih m p c' riT is the comoving density 
of electrons if the e + e~ pairs are absentj^] By using Eqn|l] 
and the condition above, we obtain 
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factor r r ~ 1, of different kinematic origin but similar 
physical properties as internal shocks, with a mechanical 
dissipation efficiency e^. These result in a proton internal 
energy, and result also in random magnetic fields with an 
efficiency e^, relativistic protons with e p , and relativistic 
electrons with ej^] In the magnetically dominated jets 
the total jet luminosity L tot in the acceleration phase 
before dissipation occurs consists of a toroidal magnetic 
field component and a proton bulk kinetic energy com- 
ponent. In the dissipation region a fraction of L to t is 
assumed to be dissipated, consuming a fraction from each 
of the toroidal field and bulk proton energy, and result- 
ing in proton internal energy and in a fraction e b which 
appears as random magnetic fields, and e p and e e which 
appear as relativistic protons and relativistic electrons. 
In both baryonic and magnetically dominated cases we 
assume eb + e p + ( e = 1, and we take ~ 0.3 and 
£b ~ 1/3 as examples in this paper. In the jet comoving 
frame the random magnetic field after the dissipation is 
parametrized by an energy density 



where 



random 



B' 2 /&tt = K eBe d L tot /{4nR 2 ph rl h c) 
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Typically, for a magnetically dominated /i = 1/3 case 
the photosphere occurs in the acceleration phase r < r sat , 

if r\ > t]t, where rjx — 150Lg3 6 r ^ & . On the other hand 
the photosphere occurs in the coasting phase r > r sat 
for rj < 77T, which is typical for baryonic cases, where 
fj, = 1 and ijT — 1900 L 53 r 7 . The Lorentz factor of 
the photosphere r p h has an r p h dependence for rj > -q T . 

being r ph = (r ph /r ) oc L^^+^-^^+^ro^ 2 ^ , 
while r p h ss 7/ in the case rj < r/x- 

In the baryonic photospheres the dissipation may be 
due to dissipation of MHD turbulence [27] or it may occur 
in the form of semi-relativistic shocks 1 15] with Lorentz 



where for semi-relativistic shocks (T IC \ ~ 1), and a 
compression ratio of k ~ 4 is assumed. For the magnet- 
ically dominated outflow, during the acceleration phase, 
the energy remaining in toroidal fields after dissipation 
is 



UL 



,toroid = (i - £ d)(i - r P h/?7)itot/(47ri?2 h r2 hC ) 



(6) 



where r p h is the bulk Lorentz factor of the protons at 
the photosphere. 

Calculations and simulations of of such baryonic and 
magnetic dissipative photospheres as well as internal 
shocks generally result in an escaping photon spectrum 
similar to the observed characteristic "Band" spectrum 
|30| . parametrized as 



1 We do not specify a particular dissipation mechanism; e.g. the 
sudden drop of photon density might trigger magnetic reconnec- 
tion as proposed by 1221 . or it could be due to MHD turbulence 
|27l or shocks 1151 . etc. We simply assume that as far as the de- 
tectable radiation the dissipation around the photosphere plays 
the largest role, and we concentrate on the photospheric dissipa- 
tion region 

2 We also assume that the dissipation gives rise to a prompt photon 
spectrum of the observed Band function type, without specify- 
ing the underlying mechanism, e.g. seed photons scattered by 
electrons associated with turbulent Alfven waves, synchrotron 
radiation from Fermi-I accelerated electron or collisional mecha- 
nism by decoupled proton and neutron | 17l .22] 1231 1271 

3 The presence of pairs will increase the radius of the photosphere 
by a factor of a few in a magnetized photosphere 24] 1281 or in a 
baryonic photosphere where dissipation is via MHD turbulence, 
or by a factor ~ 20 — 30 if baryonic dissipation is via pn colissions 
|17) . In i jv| we exemplify the effects of the effective photospheric 
radius being larger. 



dN^/dE oc (E/E hl ) x » h (7) 

in the observer frame. Observationally, for average 
bursts at redshifts z ~ 2 the mean values are £"b r ~ 300 
keV, a; p h = —1 below E^ and x p h = —2 above E^. 
In a photosphere this spectral shape is the product of 
the modification of a thermal spectrum by the dissipa- 
tion. For the purposes of this article, we treat this pho- 
ton spectrum as the input for our calculations, trans- 
formed to the rest frame of the outflow. While the bulk 
Lorentz factors in the photosphere and internal shock 



4 An alternative baryonic dissipation involves pn collisions 1171 (see 
also 1291 ) ; here for simplicity and for intercomparison with other 
models we just assume shock dissipation in the photosphere, 
whose effects are comparable to those of magnetic dissipation. 
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models may differ, for the purposes of comparison we 
adopt here as a test case the same comoving frame photon 
spectral break energy for the dissipation zones of the var- 
ious models considered, E' hl = 0.01 MeV, and x p ^ = — 1 
below EL X and x p h = — 2 above E' hr . The lower and up- 
per branches can have cut-off energies, e.g. determined 
by synchrotron self-absorption below and acceleration re- 
strictions or 77 — > e + e~ pair production above, the cut- 
off values depending on the specific model and its pa- 
rameters. For simplicity, here we adopt the same con- 
stant values of a lower limit E' min = 1 eV and an upper 
limit E' m&x = 0.5 MeV, which are adequate for our pur- 
poses since the neutrino results are insensitive to these 
values. The total luminosity of this Band-function spec- 
trum is normalized to e^Ceitotj where L tot represents the 
total luminosity. (An additional softer thermal spectral 
component can also be present at the photosphere. How- 
ever, the temperature of this component is estimated as 
T ~ o(l) X keV at the photosphere [23], corresponding to 
a thermal luminosity ^thermal ~ o(l) x 10 49 erg/s which is 
low compared to L tot and L 7 . Hence we have neglected 
this component for the purposes of the present neutrino 
calculation.) 

When the outflow encounters the external medium, it 
starts to decelerate at a radius 

D i / 31/tot^dur n1/3 

4Trni SM m p c^ri z 

= 2.54 x IO^^WIOs^Vsm^oo 73 (8) 

where an external shock forms which is also able to pro- 
duce neutrinos. Here we have assumed a uniform inter- 
stellar medium of particle density rtisM = 10 2 nisM,2crn~ 3 
and a jet outflow duration time idur in the central engine 
frame. The interstellar density value does not affect the 
photospheric or internal shock neutrinos, but it does af- 
fect the external shock neutrinos. Here we have adopted a 
density which is optimistic for the external shock neutri- 
nos, since even so the external neutrino fluxes predicted 
are low and more moderate densities such as the typically 
used nisM = lcm~ 3 would lead to even smaller external 
shock neutrino fluxes. The corresponding deceleration 
timescale is estimated as t' d ~ Rdjcq. At this radius 
deceleration Rd the external shock has fully developed, 
consisting of a forward shock, and possibly also a reverse 
shock (if the magnetization parameter a is or has become 
low enough at this radius). If present, for our parameters 
the reverse shock is marginally in the so-called thin-shell 
regime, the reverse shock having become semi-relativistic 
as it crosses the ejecta at about the deceleration time td- 

The turbulent magnetic fields generated in these exter- 
nal shocks lead to synchrotron radiation, as well as syn- 
chrotron self-Compton (SSC) and external inverse Comp- 
ton (EIC) scattering of non-thermal photons from the 
dissipation region near the photosphere or the internal 
shocks. The detailed method of calculation of these pho- 
ton spectra are discussed in [24, 3TJ[32]. As shown below, 
however, the neutrino fluence from the external shock re- 
gion is a few orders of magnitude lower than that from 
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FIG. 1. Photon spectrum in the observer frame for a typi- 
cal GRB with parameters L to t = 5 x 10 52 erg/s, ta ur = 10 
s, r] = 300 , z — 1.0, nisM = 100 cm -3 , td = 0.3, 
e c = e p = £b = 1/3, where £d is defined as the total dis- 
sipated energy from jet total energy, in the forms of e c ,e p 
and £b- The subindex "ph" refers to photosphere, "prompt" 
refers to the prompt emission from the photospheric region, 
for which a Band-like spectrum is assumed. The SYN,R 
and SYN,F are the synchrotron from the reverse and forward 
shock; SSC is the synchrotron-self Compton spectrum, EIC 
is the inverse Compton scattering of the prompt photons in 
the external shock region. (At very high energies, a Klein- 
Nishina break may be expected; however, for the neutrino 
calculation this contribution can be neglected since it con- 
tains very few photons, hence these KN breaks are not shown 
here) . A smaller dissipation value for external shock emission 
used , £b,fs = £b,rs = 0.02eb, suitable to explain the external 
shock photon emission |24| . 

the photospheric or baryonic internal shock regions, due 
to a much lower photon density leading to a lower inter- 
action rate and lower pion production efficiency. Under 
the assumptions made here, the input photon spectrum 
of a GRB with typical parameters is shown in Fig[T] as 
an exampl^] For the reverse shock, we include both the 
self-generated photons from the reverse shock and the 
prompt emission as target photons for inverse Compton 
scattering as well as for jry interaction. For the forward 
shock, we include the forward shock (FS), reverse shock 
(RS) and prompt photons. 

III. METHOD OF CALCULATION 

In the baryonic dissipation regions, whether these are 
in the photosphere or in internal shocks beyond the pho- 
tosphere, it is usually assumed that protons, as well 



5 The photospheric spectrum here does not include the effect of the 
relativistic leptons injected if we had included nuclear collisions 
1171 : the effect would be to extend the upper branch of the Band 
photon spectrum into the GeV range; however, it is the photons 
around the Band peak that affect significantly the photo-pion 
neutrino production discussed here. 
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as electrons, are accelerated through a Fermi-first order 
(Fermi-I) acceleration mechanism in the disordered mag- 
netic fields created in the region. A process similar to 
Fermi acceleration is also expected in magnetic recon- 
nection regions where layers of magnetic field of opposite 
polarity meet and drive converging flows. [33-35J. The 
particles bounce back and forth in the converging flow 
between the layers and can reach similar maximum ener- 
gies as in the usual Fermi mechanism. We assume that 
the injection of accelerated protons has a spectrum 



dN p /dE oc E~ x " (E p . 



<E n <E n 



(9) 



with x p = 2 as a nominal value. The acceleration 



timescale is t' 



r 



cycle 



iv r ' g l c = ipK/eB'c , where 



is the average gyroradius. Here we have adopted a mini- 
mum injection energy of protons E' min = 10 GeV since 
the neutrino spectrum is insensitive to this value. We 
have also assumed a high compression ratio and weakly 
disordered magnetic fields, corresponding to £ p ~ 10. 
The accelerated proton spectral energy (eqn{9| here is 
normalized to a fraction of the jet total luminosity e p L tot , 
the value of e p being discussed in section |IV| 

The maximum proton energy is constrained by the gy- 
roradius being smaller than the size of the acceleration re- 
gion r' g < i? p h/r ph , or by radiative cooling t Py&cc < i P:COO i 
where i p . CO oi is the total cooling timescale for the proton 

p,cool = *P7 + Vp + l ~BH + tsy + *7c + l ad ■ The termS 011 

the right hand side are the photohadronic, pp collisional, 
Bethe-Heitler (photopair), proton synchrotron, inverse 
Compton and adiabatic inverse cooling timescales in the 
fluid comoving frame (we have dropped the "prime" su- 
perscript here), given respectively by 



t 



■n ph (E) / dee<r p7 (e)K pl (e) (10) 

-/e TH 



t pp = cn p a pp {-f p )K pp ( lp ) (11) 
*Bif ~ 7 ( W !5 \ afCTTC f d lel - 2 n ph (-f e m e c 2 ) 

x {(2 7p 7 e ) 3/2 [log(2 7p7e ) - 2/3] + 2/3} (12) 
i- 1 = Aa T ml lp {B 2 /^)/imlc (13) 

F(E, 7p )n ph (E) (14) 

(15) 



_j 3(m e c 2 ) 2 aTC 

hc = i6 7 2( 7p -i)/3 p y w 



Q 1 « rc/fl 



where each individual cooling inverse timescale is de- 
fined as i _1 = —(d 7p /dt)j p and n p h(E) = dN/dEdV is 
the photon differential spectral density. Neutrinos result 
mainly from charged pion and kaon decays, to the first 
and second leading order of approximation respectively 
here. These charged mesons come from p 7 and pp in- 
teractions (eqn 10|11 1 . The cross section and inelasticity 
in the former channel, considering the lower threshold 
£th = 0.2 GeV, are approximated by two step-functions: 



CT p7 (e) 



340 [ib 0.2 < e/GeV < 1.0 
120 [ib e > 1.0 

'0.2 0.2<e<1.0 
0.6 e>1.0 



(16) 
(17) 



where e is the photon energy in the proton comoving 
frame. For 0.2 < e < 1.0 GeV the cross section is dom- 
inated by resonances while for e > 1.0 GeV multi-pion 
production takes over. In the single-pion resonance chan- 
nel, 7r + and 7T° are created at approximately the same 
rate. In the multi-pion channel, we assume that pions 
are created with an average multiplicity of 3, and in a 
first-order approximation the ir + , tt~ and ir° come in 
equal numbers. (For more details, see e.g. [36]). For pp 
interactions, the (thermal) protons which are not accel- 
erated in the fluid are the targets. In the fluid comoving 
frame, where the target protons are basically at rest, the 
incident proton with energy E p (E p 3> 1 GeV) has a cross 
section approximated by 



CT PP «CTpp,inci w 30[0.95 + 0.06 \og(E p )} mb (18) 
with charged pion multiplicity approximated by 



X(s)«1.17 + 0.3 logs + 0.13 log 2 s 



(19) 



in which s is the invariant energy y/s squared of the 
binary particle system. The detailed pion spectra are 
discussed in e.g. |37[l38]. However, in this paper we take 
the approximation that the pions are created at rest in 
the CM frame of the binary particle system. The error 
caused by this approximation is much reduced in the cal- 
culation of a broad spectrum in the high energy regimes 
(compared with the case in [37 or [38j),and is smaller 
than the astrophysical uncertainties in this paper. 

To derive eqn(l2] we have used a cross section rj 
(7/67r)a/(TT log(e^/2m e c 2 ) where ctf = 1/137 is the fine 
structure constant, ut — 665mb is the Thomson cross 
section and e is the photon energy in the proton rest 
frame. We note that although p 7 — >pe has a larger cross 
section than the photopion process, the effective inelas- 
ticity of the proton is smaller and the relative photo- 
pion and photopair energy loss rate for protons interact- 
ing with the peak of the vF v target photon spectrum is 
K^a^l K^a^f, ~ 100. An expression for the function 
F(E, 7 P ) in eqn 14 is given by With the above ana- 
lytical approximate expressions we can calculate the en- 
ergy fraction from the parent proton spectrum converted 
into pions, f P1 =tpi /t ? ^ cool &nd f pp =t~ p yt p ^ ooV With the 
leading order approximation that the pions are created 
at the rest frame of the protons, we can obtain the pion 
spectrum. We also rouphly approximate the produced 
Kaon number density as 1 ~ 10% of the pions from pr/ 
or pp interactions, motivated by [401 141| or simulations 
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FIG. 2. Proton inverse cooling timescales at the photosphere 
as a function of proton energy in the magnetic reconnection re- 
gion comoving frame, defined as i _1 = \Edt/dE\ in s _1 . The 
symbols used are Triangle(pg) : p-gamma (photopion) inter- 
action; Circle(bh): Bethe-Heitler (photopair); Square(pp): pp 
interaction; Dot(sy): proton synchrotron; Cross (ic) : proton 
inverse Compton; Solid line(ad): adiabatic cooling; Plus(ac) 
: Fermi acceleration timescale; Minus(TOTAL) : total cool- 
ing timescale. The astrophysical parameters for this GRB are 
itot = 10 53 ' 5 erg/s, r) — 300. At lower energies, protons are 
mainly cooled by adiabatic expansion and pp collisions. At in- 
termediate energies in this figure, protons are mainly cooled 
by photopion production and at higher energies by proton 
synchrotron radiation. 

using PYTHIA-8. Neutrinos from Kaon decays are gen- 
erally subdominant but they become the main compo- 
nent at the high end of the neutrino spectrum. At these 
energies charged Kaons suffer less from radiative cooling 
than charged Pions due to their larger mass and shorter 
lifetime. The various cooling timescales for a GRB with 
typical parameters are shown in Fig(2j 

The pion decay kinematics are well established. Neu- 
trinos result mainly from the following channel: 



being given by the condition t coo \ > idyn- The resultant 
values are summarized in Table HI 



IV. PARAMETERS AND NEUTRINOS FROM A 
SINGLE SOURCE 

We numerically compute the neutrino spectrum by us- 
ing the method described in the previous section. Several 
parameters are needed: the jet total luminosity L to t and 
its duration in the source frame tdur, the target photon 
luminosity L 7 = e^e-^totj the Fermi accelerated proton 
luminosity and its power-law spectral index x p , the mag- 
netic field calculated from the parameter €b and eqn|5] 
the dissipation radius R, the outflow energy to mass ratio 
•q and finally the source redshift z. 

We discuss three main representative scenarios: an ex- 
treme magnetic photosphere model where V oc r 1 / 3 , a 
baryonic photosphere model where T oc r, and a modi- 
fied internal shock (IS) scenario where two ejecta shells 
of different bulk Lorentz factors collide. For the modified 
internal shocks we assume that a high mechanical dissi- 
pation efficiency is achieved, e.g. [5lUl|, and that the 
photon spectral issues raised about traditional internal 
shocks are avoided in such mechanisms. As far as neu- 
trino production, the location and seed photon spectrum 
is similar to that in the standard internal shock, but with- 
out the simplifications of H] in the neutrino physics; 
that is, we treat the pj interaction with the whole pho- 
ton spectrum, not just the break region, and include be- 
sides the A-resonance also multi-pion effects, Kaons and 
detailed secondary particle distributions for the charged 
meson and muon decay, as in e.g. |431 144| . We assume 
for the internal shocks a dissipation efficiency = 0.3, 
which for comparison is taken to be similar to that of 
the photospheric models. For the two photospheric mod- 
els, the extreme magnetic one satisfies rj > r)x, while for 
the baryonic one r\ < r/x, so the photospheric dissipation 
radii are 



tt* -> /i* + v^Vy) ->■ e ± + v e (v e ) + + v u 



(20) 



which we calculate in detail following the method in 
|42| . A high energy pion may lose a significant fraction 
of its energy through synchrotron radiation before it de- 
cays. Therefore we first calculate the pion spectrum af- 
ter the cooling process has set in (using a method similar 
to eqns. (l3|l4] and[l5]). Then we calculate the neutrino 
spectrum from the 'final' pion and muon spectra. The 
fj, has a longer mean life-time and smaller mass which 
makes its synchrotron cooling more severe than that of 
charged pions. Finally we note that the leading decay 
channel of the charged kaon is the same as that of the 
charged pion so in this sense they can be viewed as "ef- 
fective pions". 

While the maximum proton energy at injection is de- 
termined by the condition t acc < i CO oi; the maximum 
cosmic ray energy of the escaped protons can be smaller, 



*ph _/4 8/ V 3/5 



''o 



for /i = 1/3; 
for fi = 1 , 



(21) 



where t\t is defined in eq. Q. For the internal shock 
models, the dissipation radius at which these shocks oc- 



cur is 



CT] ^var ■ 



(22) 



Here r\ is the average Lorentz factor of the two shells 
and £ V ar ^ 1 nis is the variability timescale which rep- 
resent the time interval between the ejection of the two 
shells in the source frame. Note that a range of t var 
is indicated by observations, extending down to i var = 
0.001 ~ 1 s. This introduces a large uncertainty in the 
internal shock radius and the corresponding final neu- 
trino spectrum. Here we use optimistic values for the 
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TABLE I. Parameter list for different models calculated. The common parameters are: jet total luminosity L to t = 10 
erg/s, source frame duration t^ur = 10 s. source redshift z = 0.1, and the dissipation partition fractions e c = eb = e P = 1/3 
with ed = 0.3 . The first column identifies the different curve symbols and the figure used. The M and B in the second column 
refer to "Magnetic dominated" and "Baryonic dominated". The third column identifies the type of dissipation reg ion : Ph for 
photosphere with i? p h from eqn|3l 



= Ph for pair-photosphere, and IS for internal shock of radius Ris from eqn. 
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next columns R = 10 Ris cm gives the corresponding radii, 77100 is the initial dimensionless entropy, TiooQi? and B 5 



In the 

B 5 @R 



are the bulk Lorentz factor and comoving magnetic field (in the unit of 10 G) in the dissipation region. The max -£/CR,PeV is 
the maximum escaping cosmic ray (proton) energy in the source frame, calculated by setting t pl = td yn , or t p , sy n = £d yn , or 
i P ,cooi = £ P ,acc, whichever gives the smallest E. The number of /^-events are in the last two columns; these are estimated from 
the neutrino flux and the effective area of the IceCube 86-string configuration, for a source at z — 0.1 or z = 1.0. 



internal shock neutrino production, for comparison pur- 
poses. Unless specified otherwise, in the following we 
assume a nominal parameter set of L tot = 10 53 5 erg/s , 
^dur = 10 s (source frame) , r\ = 300 , e^e = 0.1 cor- 
responding to an isotropic equivalent total photon lumi- 
nosity of L 7 w lO 52 5 erg/s (which is roughly the average 
luminosity from GRB statistics). The exact energy parti- 
tion fractions in the jet are not well known, here we have 
assumed = e e = e p = 0.33, td = 0.3. As an example 
of the fluence in the observer frame , we consider a GRB 
at redshift z = 1 corresponding to a luminosity distance 
of 6.6 Gpc, and at z = 0.1 corresponding to a luminos- 
ity distance of 450 Mpc in a standard ACDM cosmology 
with 

fl m = 0.28, fl v = 0.72, H = 72km/s/Mpc. (23) 

The detailed parameters of different models are listed in 
Tablelllfor which the neutrino spectra are plotted in figs. 
Sand M 

The neutrino spectral fiuence (iV+Sv) from a magnet- 
ically dominated fi — 1/3 GRB is shown in Figj3] includ- 
ing both photospheric and external shock contributions. 
In principle, the neutrino flavor distribution at the source 
can be calculated from by Eqn |20| and then recomputed 
at the observer frame after neutrino oscillations. However 
the large range of distances, source geometry and density 
distribution introduce a large variability in final result, 
so here simply approximate the received neutrino flux as 
having equal numbers in all three flavors. In this type 



of models the dominant neutrino emission comes from 
the magnetic photosphere. The neutrino spectrum from 
the external shock peaks at a higher energy because the 
magnetic fields and photon densities there are lower than 
in the photosphere, and the charged mesons suffer much 
less synchrotron cooling before they decay. A uniformly 
distributed interstellar density of njgM = 100 cm -3 is as- 
sumed for the external shock calculation. The neutrino 
fiuence is nonetheless low compared to the photospheric 
fiuence, due to the low target photon and proton column 
density and therefore the low pp or pry collision rate. A 
smaller value of nisM would lead to a more distant shock, 
and even lower neutrino fluences. 

Fig(3j as well as Fig |4|5| also illustrate the effect on 
the neutrino spectra of including additional physical pro- 
cesses besides the pj production from the A + -resonancc 
used in many previous studies, including the recent Ice- 
Cube GRB data analyses |J. The processes included 
in Figs{3] [4] and the rest are the production of ^ via 
A + in as well as pj, as well as multi-pion production 
and K ± in pp and pj. We note that pp collisions can 
become important at lower energies (before pj interac- 
tions set in). This is especially true if the dissipation 
radius i?dis is small, where the pp collision optical depth 
T pp approaches unity or above. We also note that the 
neutrino fiuence is comparable to the photon fiuence in 
fig^ This may be a relatively conservative value; in 
some other works a higher acceleration efficiency , pio- 
nization efficiency and proton luminosity e.g L p = 10 L 7 
are adopted corresponding to a neutrino luminosity L v ~ 
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FIG. 3. The observer frame ^ M + ^ M fiuence E 2 dN/dE for 
a jit = 1/3 magnetically dominated burst at z = 1. Solid 
line(i^ ph ): total muon neutrino spectrum from A-resonance 
(dotted), P7 multi-pion production (circle), pp collision (cross) 
and Kaons (square). Large dashed line (v RS and v FS ): total 
neutrino spectrum from external reverse and forward shock 
in the early afterglow. For the forward shock neutrino cal- 
culation, we included the photons from the reverse shock as 
well. 



FIG. 5. The neutrino spectrum from a magnetic photosphere 
(red upstanding triangle), a baryonic photosphere (blue in- 
verted triangle) and an internal shock model (black square), 
for T) = 300 (filled) n = 1000 (empty). Both the two lines 
for internal shock models have used t vaT — 1 ms. The other 
parameters are the same as in fig|4] Note that a slower jet 
produces a higher neutrino flux in all these cases. The more 
detailed quantities for those models (and in fig|4| are listed 
in Table HI 
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FIG. 4. The + fiuence spectra for different models 
assuming a source at redshift z = 1. Ltot = 10 53 ' 5 erg/s, 
r\ — 300 are used. (External shock neutrino spectrum is not 
calculated here.) The models shown here are the fi — 1/3 
magnetic photosphere (M,Ph), fj, = 1 baryonic photosphere 
(B,Ph) , pair-photosphere dissipation (M,e ± Ph or M,e ± Ph) 
and internal shock shock models (IS) with f var = 1, 10, 100 ms 
(black, dark gray, light gray) 



a few L 7 . 

In Fig(4] we also show the possible effect of having a 
dissipative photosphere at a radius larger than that given 
by Eqn.([3|, e.g. due to pair effects. Without specifiyng 
an underlying mechanism giving rise to the dissipation 
and Band function spectrum this increase is uncertain. 
If the photon spectrum is generated by the scattering of 
thermal electrons associated with turbulent Alfven waves 
|27| . the maximum energy of the comoving photons can 
hardly exceed the electron rest mass m e c 2 and negligi- 
ble pairs are produced. On the other hand if the prompt 
emission is due to synchrotron radiation from Fermi-I ac- 



celerated electrons e.g. [35] and the attenuation beyond 
E' max is due to 77 —> e + e~ process, we can calculate r 77 
and the amount of pairs. The existence of pairs with such 
parameters (Fig[T]) could boost the pair-photosphere ra- 
dius to ~ 20 times the original i? p h. In the collisional 
scenario where protons and neutrons decouple |17| . high 
energy 7s and injected poisitrons from decaying pions in- 
duce leptonic cascades and significant pair production, 
which gives a similar boost of ~ 20 to a pair-photosphere 
radius. Without going into specific model details, we just 
consider for simplicity the dissipative and spectrum for- 
mation effects associated with such larger effective pair 
photospheres, and use this radius for calculating the neu- 
trino spectrum. In Fig|4]the curve labelled M,e ± Ph is for 
a magnetized dynamics photosphere a factor ~ 20 larger, 
and the curve labelled B,e ± Ph is for a baryonic photo- 
sphere a factor ^3.6 larger than the value of Eqn.(|3]). 

The case of an internal shock is also shown, for a dis- 
sipation radius estimated by Eqn |22| and variability time 
^var = 1 ms. Qualitatively, the harder spectra from larger 
radius dissipation regions arise because the magnetic field 
and photon density is lower at larger radii, reducing the 
pion and muon electromagnetic cooling. 

The effect of a magnetic versus a baryonic jet bulk dy- 
namics on the neutrino signatures are compared in Fig|5] 
where we plot the neutrino spectra from a magnetic pho- 
tosphere where T oc r 1 / 3 , a baryonic photosphere where 
r oc r and an internal shock model. Each case is calcu- 
lated for two different terminal Lorentz factors, i] = 300 
and r) = 1000. 

The neutrino fiuence in general decreases for increas- 
ing 77 in the baryonic photosphere model. The analytical 
expression for the photospheric radius is R cx rj~ 3 in this 
model (Eqn(3]), and having assumed a fixed Band func- 
tion in the comoving frame, the number density of target 
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photons is n' ph ~ u' vh /E'AE' cx L/R 2 rf cx rf . There 



ph 



fore, the inverse cooling timescale t^" 1 oc n' ph cx 77 4 
We also have i^~ 4 ~ ?7c/i? ~ rf ~ ^ . The pionization 
efficiency / ff = t'^ 1 /t'~^ has no dependence on 77 in the 
leading order approximation, or on R, which is treated 
as a function of 77 here. But, the synchrotron cooling for 
n ± , fi^ (t aync /t dcc cx 77~ 4 for given tt,/j energy) dimin- 
ishes the neutrino flux for larger 77. From Fig [5] we see 
that the peak energy of 77 = 1000 is smaller than 77 = 300 
case in the comoving frame. This indicates strong pion 
cooling taking place. Therefore a higher 77 is associated 
with a smaller v flux. 

For the magnetic photosphere, a larger 77 gives smaller 
.Rph and T p h. Both are advantageous for neutrino pro- 
duction (although a higher magnetic field lead to stronger 
cooling of pions and muons, it is less severe than the bary- 
onic case, e.g. see B$@R column in Table]!}. 

In Figj6] upper panel, we compare the luminosity de- 
pendence of the single source spectra for a standard 
77 = 300 of an "optimistic" internal shock model (with 
a high dissipation rate and i var = 1, 10, 100 ms), a 
magnetic photosphere model and a baryonic photosphere 
model. The optimistic internal shock model with i var = 1 
ms gives the highest flux at energies < 100 TeV, while the 
magnetic photosphere models have a lower flux at these 
energies, with their spectrum also peaking towards higher 
energies. In Fig(7] upper panel we show those cases with 
7/ = 1000. We see that the magnetic photosphere case is 
the least affected by 77. 

An interesting possibility in the case of magnetic dis- 
sipation regions, e.g. due to reconnection, is that they 
may produce an an accelerated proton spectrum which 
is harder than the typical Fermi case of dN/dE cx E~ 2 . 
As discussed by |4"5l HB] for an extreme case all the pro- 
tons entering the acceleration process are essentially con- 
fined by the magnetic field and the zero escape probabil- 
ity leads to a spectrum dN/dE cx E^ 1 . If this scenario 
is valid, the bulk of the energy for the accelerated pro- 
tons is concentrated in the high energy end of their in- 
jection spectrum. The neutrino energy associated with 
these protons lies in the PeV-EeV energy range. This en- 
ergy is in the sensitivity range of IceCube and proposed 
AM ANN A neutrino detector. However this scenario re- 
quires significant magnetic reconnection process where 
the toroidal magnetic field is still present. Charged pions 
and muons at this energy suffer strong synchrotron cool- 
ing which suppresses the neutrino spectrum significantly. 
Therefore, no significant neutrino emission is expected 
from this scenario (since at lower energies there are too 
few protons due to the nature of a dN/dE cx E^ 1 spec- 
trum.) 



V. DIFFUSE NEUTRINO BACKGROUND 
FROM VARIOUS GRB MODELS AND 
IMPLICATIONS 

Since the individual source fluxes are very low, except 
for the unlikely event of an extremely nearby occurrence, 
it is useful to consider the cumulative diffuse neutrino 
flux from all GRBs in the sky. We calculate this diffuse 
flux based on two methods. Method-I uses a GRB lumi- 
nosity distribution (luminosity function) and a redshift 
distribution |47| . given, respectively, by 



< L 7 < L» 



(L y /L*y 



Rgrb(z) cx 



(i + zY 

(1 + 2)' 



Z < Z\ 

z > z\ 



(24) 



(25) 



Here L 7 is the peak photon luminosity (here mostly in the 
0.1-1 MeV range) and we have used the following values: 
L min = 10 50 erg/s, L. = 10 52 5 erg/s , L max = 10 54 erg/s, 
mi = -0.17 , m 2 = -1.44 , Tii = 2.07 , ti 2 = -1.36 , 
Z\ = 3.1 which are best fit values used in [47 . The 
differential comoving rate of GRBs at a redshift z is 



R(z) 



R, 



GRB 



(z) dV 



(1 + z) dz 



(26) 



where V(z) is the comoving volume in the standard 
ACDM cosmology given by parameters in Eqn|23] and 
the factor (1+z) accounts for time dilation effect due to 
cosmic expansion. The differential number of GRB per 
unit redshift is given by 



&N(L 1 ,z) = pQ4>{L 1 )R{z)A\ogL 1 dz 



(27) 



where po = 1.3/yr/Gpc 3 is the local rat^jof high lumi- 
nosity GRBs (not including a separate family of objects 
called low luminosity GRBs (48]). 

An alternative method (Method-II) is to use the Wis- 
consin GRB catalog |http: / / grbweb.icecube.wisc.edu/ 



We consider only long GRBs and use the parameters 
such as the Band photon index, fluence, redshift, pho- 
ton peak energy etc. to compute the neutrino spectrum 
from each individual burst and stack them together. Fi- 
nally we normalize the resultant flux to an all sky rate of 
700 GRBs/yi[^] These two methods are then used to cal- 
culate the diffuse neutrino flux assuming that the GRB 



the proton radiative cooling is also relevant; however, it is rela- 
tively weak due to their heavy mass compared to pion and muons 



7 if we use the simplest detection criteria for SWIFT dNp^/dt > 
0.4/cm 2 /s, the number of GRBs which meet this criteria is about 
600/ yr which is a reasonable all sky rate. 

8 Since these instruments have a limited sky coverage and opera- 
tion time bin, a smaller number of GRBs per year are actually 
recorded 
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FIG. 6. Upper panel: Neutrino fluence from a single GRB 
from different dissipation regions. Red, dashed: magnetic 
photosphere; blue, dotted: baryonic photosphere; Dot-dash: 
baryonic internal shock for t var = 1 , 10, 100 ms. (black, 
dark gray, light gray). These are computed for different 
luminosities (in each model, L 7 = 10 53 , 10 52 , 10 51 erg/s 
(top, middle, bottom). Lower panel: Diffuse + P p neutrino 
spectral flux from the three models above (same line style), 
calculated for an all-sky GRB rate of 700/yr using statis- 
tical Method I (thick lines) and Method II (thin lines; see 
jjyj. Also shown is the IceCube collaboration's representation 
of the diffuse flux from a standard Waxman-Bahcall internal 
shock model, and the IC 40+59 observational upper limit (see 
Fig. 3 of [4 for description). The gray zone labeled ATM is the 
atmospheric neutrino spectrum. The plots from both panels 
suggest that the occasional electromagnetically bright GRBs 
can contribute significantly to the total diffuse flux. 



neutrinos are due to a magnetic photosphere, a bary onic 
photosphere, and an internal shock model (Fig |6|7[ ). In 
Fig(6] we plot those cases with rj = 300 while in FigjT] 
77 = 1000. 

An inspection of this Fig 6]7 (lower panel) shows that 
the above models, with the parameters used in this paper, 
predict a diffuse neutrino flux which is likely to be within 
the current constraints set by IceCube 40+59 string ob- 
servations, as suggested by the fact that they lie below 
the two IceCube constraint lines labeled WB for IC and 
IC 40+59. A caveat is that these constraints are upper 
limits on the flux as a function of the break energy assum- 
ing a Band function with a = —1,0 = —2, or assuming a 
slope -2. They are not directly applicable to other types 



FIG. 7. The parameters and conventions are the same as 
those in fig[6] except we have used r\ = 1000 (in figjfjjry = 300 
is used). 



of spectra; thus, these constraint lines are here intended 
only for a rough comparison; they would have to be re- 
evaluated for the spectra shown here. Nonetheless, they 
do provide some guidance, and it appears that the only 
models which are close to being constrained at present 
are those using the most optimistic parameters; e.g. for 
internal shock models with a i V ar = 1 ms one would ex- 
pect a small dissipation radius leading to a high r pp and 
a high neutrino fluence in the lower energy range from 
pp collisions (where r p7 has saturated to unity with t pp 
increasing, making pp collisions more important.) Also, 
if we were to assume a higher value of e p /e e , such as 
10 (in this paper we used a value of 1, see also section. 
IV) and/or if we were to adopt a lower magnetic field 
fraction (eg ~ 10~ 2 , rather than the e^B = 0.1 used 
here), the diffuse neutrino fluence would be likely to vi- 
olate the above constraints, especially for the internal 
shock model. However, this is for the optimistic internal 



shock case where one uses t v 



1 ms. whereas there is 



a larger uncertainty in this quantity, and in most bursts 
i V ar can often be several orders or magnitude larger (see 
comments in the third paragraph of \ IV I. For such larger 
(and more reasonable) values of i var the model would ap- 
pear to be still compatible with the constraints, even if 
a high e p /e e ratio were assumed. 

We have also included the approximate atmospheric 
neutrino background in Fig|6|7| It is worth noting that 
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GRBs are transient sources whose prompt emission du- 
ration in the observer frame is within tf^. ~ 100 s. The 
angular resolution for TeV neutrino and above is within 
1 deg 2 . Therefore, considering the search time bin and 
the small solid angle set by optical observations, the ef- 
fective atmospheric neutrino background is well below 
these GRB diffuse fluxes. In other words, even one or 
two muon events in IceCube correlated with photon de- 
tections would give a high signal to noise ratio. 

A caveat for such calculations of the diffuse back- 
ground is that the statistical description of both the 
source physics parameters and the spatial-temporal dis- 
tribution of the sources has large uncertainties; this ap- 
plies to both the luminosity function (Method I) and the 
observed burst catalog (Method II). In Mcthod-I, the pa- 
rameters in eqns. (24] 25 1 have large uncertainties; es- 
pecially in the low redshift region (e.g. z < 0.3), very 
few GRBs are observed. However, in order to expect 
more than one muon event in IceCube, we need a GRB 
of moderate or high luminosity located at low redshift 
(such as the example GRB in figHwith L 1 = 10 52 5 erg/s 
and z — 0.1, which results in about 0.2 muon events in 
the IceCube 86-string configuration.) The estimated rate 
of GRBs which satisfy iV M! i cecu b e > 1 is about o(l) per 
ten years. In Method-II, we actually have a very lim- 
ited number of GRBs with well measured redshifts in 
the catalog. Those without a redshift are assigned a de- 
fault redshift value of z = 2.0. Thus, improvements in 
the photon-based statistics of bursts, as well as neutrino 
observations over multi-year periods appear required. 



VI. DISCUSSION 

We have calculated the + signatures expected 
from photosphcric GRB models where the dynamics is 
either magnetically or baryonically dominated, including 
also the effects of external shocks, and have compared 
these signatures with those expected from the baryonic 
internal shock models. This comparison is timely in view 
of recent developments, the most pressing being the re- 
cently published IceCube constraints |3J 0] on the "stan- 
dard" internal shock GRB models. The exploration of 
alternative models to the internal shocks has its own 
separate motivation, independently of the IceCube ob- 
servations. One reason is the increased realization that 
magnetic fields may play a dominant role in the GRB 
phenomenon, and the dynamics of magnetically dom- 
inated GRB jet models |^] are being considered in de- 
tail [2D1 [50] • Also, issues related to the efficiency 
and spectrum of standard internal shock models of the 
prompt 7-ray emission [BHE] have led, on the one hand, to 
considering modified internal shock models that address 



these issues, e.g. [101 [TTl [ST] . etc., and, on the other hand, 
to considering the photospheric emission as the source of 
the prompt 7-rays, including both baryonic photospheres 
|13ffT5] and magnetic photospheres |23l |2"4"1 152"] . Some 
previous calculations of neutrino spectra from baryon- 
ically dominated photospheres have been carried out 
[531 154"] , but so far none from dissipative photosphere^"] 
or magnetic photospheres, which we treat here 

One of the notable features about the neutrino emis- 
sion from dissipative photosphere models, both in mag- 
netized of baryonic dynamics, is that the peak energy 
of the spectrum (which is around 0.5-1 PeV, as seen 
in Fig. [3]) differs from the peak energies of a typical 
Waxman-Bahcall (WB97)[T] simplified standard internal 
shock model such as used in the recent IceCube studies 
(and in j5] I56|). The peak energy is also higher than 
those for the internal shock model with £ var = 1ms and 
7] = 300, but lower than those with higher i var and r\ val- 
ues. In part, this is due to the inclusion, in addition to 
A resonance production, also of multi-pion effects, Kaon 
decay, etc., which also gives a steeper spectrum above the 
peak, compared to the flat slope of E 2 dN/dE above the 
peak expected in the usual WB97 internal shock model. 
This steepening of the spectrum above the peak is not 
unique to the photospheric dissipation model, it occurs 
also in internal shocks when we include the above addi- 
tional physics beyond the A-resonance, as can be seen 
in Fig. [1] (Such a steepening for internal shocks is also 
found bv p]V 

The energy at which the spectral peak appears de- 
pends generically on the radius of the dissipative and 
photon escape zones, which here we have assumed to be 
colocated, whether it be a photosphere or an internal 
shock. The effect of different radii is illustrated in Fig. 
[4] which compares the spectra and peak energies of differ- 
ent photospheric dissipation zones at different radii and 
also two internal shocks at two different radii (or vari- 
ability times). The harder spectra from larger radii are 
expected because of the decrease of the magnetic fields 
and and photon densities, which allows the higher energy 
secondary pions and muons to decay before significant 
electromagnetic cooling has taken place. 

The difference between a magnetically dominated pho- 
tosphere (M,Ph), a baryonically dominated photosphere 
(B,Ph) and internal shock (IS) models are compared in 
Fig ]4|5| Fig ]6|7[ upper panel) and Table]!] The dissipa- 
tion region for the 'M,Ph' typically lies in the accelera- 
tion phase of the jet where there is a smaller T p h than 
'B,Ph' and 'IS' model. The neutrino spectrum in the 
observer frame is affected by the efficiency of pp and pr/ 
interaction, the cooling of the secondary charged par- 
ticles (which is mainly synchrotron, determined by the 
comoving magnetic field) and finally Lorentz boost and 
redshift. The differences in neutrino spectra from these 



Magnetically dominated jets are naturally expected if the black 
hole energy is extracted via a Blandford-Znajek mechanism |49| 
, or if the source is a temporary magnetar 1211 . 



As we were ready to submit we received a preprint on this subject 
I55| . with results compatible with ours. 
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models are not quite significant for rj — 300 where the 
dissipation region, magnetic field and Lorentz factor are 
similar. For 77 = 1000, the 'B,Ph' and 'IS' have harder 
spectra than the 'M,Ph' model. The 'B,Ph' radius here 
is rather small where pp, pj and coolings are all efficient. 
The large Lorentz boost factor finally pushes the peak 
neutrino energy in the observer frame over the 'M,Ph' 
model. For 'IS' model with rj = 1000 the dissipation 
radii is much larger than 'M,Ph' and 'B,Ph' case. The 
inefficient cooling results in the highest peak energy of 
the three models. However, the flux is also the lowest 
because pp and p'y are the least efficient of these models. 

The maximum source-frame energy of the accelerated 
protons which are able to escape the acceleration region, 
for the models considered here, are shown in the third 
from the last column of Table|TJ Leaving out any consid- 
eration of the diffuse flux, it is seen that only the inter- 
nal shock models approach the highest energies associ- 
ated with the GZK limit, as originally suggested by [lj, 
while all photospheric models fall several orders of mag- 
nitude below this energy. We have not done an exhaus- 
tive parameter searcrPM since our emphasis has been the 
neutrino emission, but it is apparent that photospheric 
models would not be competitive GZK sources. 

The detection of neutrinos from individual single 
sources with the 86 string IceCube, as seen from Table|T] 
in the last two columns, is rather difficult, the number 
of expected muon events at best being ~ 0.2 for a mag- 
netic photosphere model of average luminosity at z = 0.1. 
The only hope may be the statistically rare observation 
of a nearby bright GRB , or through observations of the 
diffuse flux. 

The diffuse flux offers higher prospects for an eventual 
detection, and a comparison of the diffuse flux expected 
from the three different types of models discussed, Figs. 
[6] and [7] lower panels, shows that an internal shock sce- 
nario with optimistic parameters such as i var = 1ms and 
77 = 300 comes closest to being constrained by the current 
IceCube limits [1]. This confirms the recent calculations 
of [431 l6"U| 1ST] . On the other hand, a magnetic photo- 
sphere model is far from being constrained by the current 
limits (although the limits will have to be re-evaluated for 
the specific spectral shapes) . A baryonic dissipative pho- 
tosphere model has an even lower flux, and would take 
the longest to be detected. Figj7] shows that a higher 
Lorentz factor, in this case r\ — 1000, decreases the con- 
straints. The diffuse flux from a magnetic photosphere 
is the least affected, and both baryonic photosphere and 
internal shock models have higher peak energy but lower 



flux. 

In conclusion, the calculations discussed in the previ- 
ous sections indicate that the current IceCube limits are 
not yet sufficient to distinguish between the dissipative 
photosheric models based on magnetic or baryonic dy- 
namics, although they are approaching meaningful lim- 
its for baryonic internal shocks outside the photospheres. 
Also, as of now it is not yet possible to set stringent 
limits on the value of e p /e e , i.e. the putative ration of 
accelerated protons to accelerated electrons. One of the 
reasons, for the photospheric models, is that the mag- 
netic reconnection or dissipation scenario is less straight- 
forward than the simple internal shock scenario, having 
both more physical model uncertainties and parameters 
associated with them, e.g. the geometry of the plasma 
layers or the striped field structures, resistivity, insta- 
bilities, reconnection rate, etc. Both for photospheric 
and internal shock models, the fraction of protons that 
are injected into the reconnection or Fermi acceleration 
process is uncertain. The acceleration timescale and es- 
cape probability depend on the geometry and dynamics 
of magnetic reconnection or acceleration region. These 
factors are also crucial for determining the injection of a 
cosmic ray proton spectrum and its effect on the neutrino 
flucnce, maximum neutrino energy and spectrum. 

For these reasons, and in preparation for future im- 
proved limits from longer observation times and from 
possible future higher sensitivity detectors, we have here 
investigated the typical neutrino spectral features and 
flux levels expected from three of the basic GRB mod- 
els currently being considered. There are reasonable 
prospects that the detection or non-detection of these 
neutrino fluxes in the next decade with IceCube or next 
generation of large neutrino detectors will shed light on 
the underlying physics of the GRB jets and on the cosmic 
ray acceleration process in then. 
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